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1. INTRODUCTION 

Various problems of Mathematical Science are often resolved within reformulation of different 
mathematical problems like integral equations. So, the study of integral equations and the ways of finding 
methods for it are quite helpful. Now-a-days, Volterra’s equations have been implemented in the fields 
of Mathematics and Technology. Some valid methods, for solving Volterra’s equations have been developed 
by many researchers. The chebyshev polynomials is one of the methods to find the solution of integral 
equation by Rahman [1]. In particular, Huang [2] has solved the unknown function by using the Taylor’s 
expansion. Mixed quadrature has faster convergence due to higher degree of precision for evaluation of real 
definite integrals. Piessen [3], Wazwaz [4], Jerri [5], Kreyszig [6], Pundir [7], Swarup [8], Conte 
and Boor [9] developed a method for computing integral transforms by Chebyshev polynomial 
approximations. Mixed quadrature has faster convergence due to higher degree of precision for evaluation of 
real definite integrals. The mixed quadrature of higher precision has been established with linear convex 
combination of Gaussian and Newtonian type rule of lower precision for single variables. Recently mixed 
quadrature has been successfully used for the numerical solution of integral equations as well as finite 
element methods. A step have been carried out by Dash and Das [10], Jena and Nayak [11] Tripathy 
et al. [12]. Das and Pradhan [13] Jena and Dash [14], Rokhlin and Wandzura [15], Jena and Nayak [16], 
Islam et al. [17], Rathod et.al. [18], Chakrabarti and Martha [19], for numerical solution of non linear 
Fredholm equation of second kind. Oliver [20], Davis and Rabinowitz [21], Dash and Das [22] proposed 
identification of some Clenshaw-Curtis rule with Fejer rules. Motivated by the excellent performance 
of these methods, we have applied the resolvent kernel to find the solution of Volterra’s integral equation 
for numerical integration of real definite integrals for single variables which has been applied in the field 
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of RLC circuit. The proposed work to find instantaneous current in RLC circuit with the help of Volterra’s 
equations with mixed rule for the integral equation of type: 


u(p)= f(p)+ A] Mp—t)ule)at, (ly 


where p € [0,t] f(p) is the source function and k is the kernel which are given and u(p) is the unknown 


function and A is anon zero real or complex parameter. 

Application of Kirchhoffs Voltage law transfoms the differential equation of RLC circuit 
to Volterra’s integral equation of second kind then it can be made higher precision rule for higher rate 
of convergence towards the exact value of the integral by mixed quadrature rule of Clenshaw Cutis seven 
point rule and Lobatto five point rule by Singh and Dash [23], Laurie [24] and Atkinson [25]. This proposed 
work is organized as follows. Section 1 meants for the introduction. In Section 2, Kirchhoff’s Voltage law 
and an approximated solution for the Volterra integral equation of second kind is established. In Section 3, 
a mixed quadrature rule is constructed which is approximately convergence to the exact solution and error 
bound is estimated for the suggested method. In Section 4 approximate solutions for various problems 
are found. In Section 5 some conclusions are drawn. 


2. BASIC PROPERTIES OF KIRCHHOFF’S VOLTAGE LAW AND VOLTERRA 
INTEGRAL EQUATION 
In this section, some basic properties have been used for definitions of Kirchhoff’s Voltage law 
and Volterra integral equation of second kind of (1). According to KVL, the integro-differential equation 
for RLC circuit is given in Figure 1. 


1 
LI +RI+— | 1dt=E(t) 2) 


where E(t) is the electromotive force. Ey =constant, I =Inductance, R =Resistance, C =Capacitance 


and 1(t)=The instantaneous current. 


C 
. : 3 L 
Ea) 
Figure 1. RLC circuit 
The integral (1) after differentiation becomes: 
ry Rpy 1 2 (3) 


L CL L 


The (3) is anon homogeneous second order differential equation. 


Bulletin of Electr Eng & Inf, Vol. 9, No. 2, April 2020 : 801 — 807 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 im) 803 


2.1. Volterra’s integral equation 

Here some basic theorems are given below to solve the Volterra’s integral equation of second kind 
which is of the form (1). 
Theorem-2. 1 





The initial value problem which is of the 


2 
os ef ane a Dt) =f (t) with the initial condition 


I(c)=1,, I'(c)= [/, is transformed to Volterra’s integral equation of second kind is of the form 
t 
1, +[1,+1C(c)] (t-c)+ [(t-p) f(p)dp- 


flee) +(-p){0(»)-c'(»)} ev) ap 


a 


I(t)= 


Proof: 





Here we have given the differential equation 


2 
+O e+ D(t)I =f (r) with the initial 
condition I(c)=I,, Ic =I. Then integrating twice with respect to f from C tot we get 


I, +[ 1) +1)C(c) | (tc) J+J p)dp- 
I (0 )= t which is the required Volterra’s Integral 
fle) +@-p){P(v)-C(v)} |p) ap 


equation of second kind. 
Theorem-2.2 


The resolvent kernel of the kernel k(p, tis a polynomial of degree (n—-1) in f which is express 


in the form  k(p,t)=D,(t)+D(t\t—p)+D,(t\t—p) +... D, Mt) (¢—p)'' +... 


(n—1) 

















1 n 
where the coefficient D, (r) are continuous in the interval (0, c| is R(p, t; A)= “7 . : o(t, D3 A) . 
ID 
Proof: 
onsider R p,t;A ee ue Q t, pA 
¢ hap (4) 
The auxiliary function ~ satisfies the following conditions: 
flO, MED iat Sie coy 6) 
dp dp" dp n-l 
therefore, the functional relation reduces to: 
d" f d" 
® — ak(p,t)+A[ kz, p)Alt,z:d)ds s 
dp : ds 


integrating by parts in (4) and (5) 
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p= £8 +f DE v()L+..0 


dp n dp n-1 





II 
i) 











hence the resolvent kernel is of the form. 


n 


R(p,:4)=-5 Wt p:) 


In this case study, we have used a polynomial in the place of the unknown function I(p) 


for the transcendental function for making: 








a) » for i=1,2 (Mrigonometric function) 
er re (7) 
iD sairoe ee 


Ft (t), if t=2 (Exponential function) 
Ft) ) if t=3 


> 


exponential function (oscillatory function) to convert smooth function. 


3. CONSTRUCTION OF MIXED QUADRATURE RULE 
In this section, a mixed quadrature rule has been constructed by the convex combination 


of Clenshaw-Curtis seven point rule Roa (f) and Gauss Lobatto-five point rule R,<(f )and also the error 


bound estimated which has been given by the following theorems. 
Theoren-3.1 


The quadrature rule and error for the smooth function f(p) which is defined on OS Dp <1 
is obtained by convex combination of Clenshaw-Curtis seven point rule Ralf) and the Gauss 


Lobatto-five point rule R,(f) rules formed mixed quadrature mule of degree of precision nine 


Recus(f), then 


22) 28a (f)+7R,5( he all 28a (f)+TE,s(f)]| where 


1 


Recah )= 13 ail Roof + TR, (f )I and Eee71s (f == [Beg (f + TE,s (f )I 


Proof: 


ea oie+ri}+s9f (2) +12} 


= 14a s(- 5) + (5) } +1640 
rdsrafrara greases {-P} [Ff ssr0 
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Each of the rule Ror (f) and R,;(f) is of degree of precision seven where Bae (f) is the error due to 


Rec lf) and E,(f) is the error due to R,(f) 
Uf)=Recf)+ Ecol) ® 
(f)=R,5(f)+ Es(f) (9) 


Expanding in (8) and (9) with the help of Taylor’s series, we get: 


ee viii 1 : 177 i 
Ecalf) = ire qcai )+ so8a07 © aor20aa7 (0)... 














a ae 128. 8256 oxi 
Elf)= 220587 (0) 377x107 (0) 15606512" (0)... 


Multiplying (2) with (8) and (3) with (9) and adding them respectively. 


Is [rarest 51 28Ree(N+7 Ris Nt 


sell28 Fecal f) +7E,(f)] 

Then [ (f = Reee )+Eceus(f ) 

where: 
Rous )= 1 Frill 2RReal )+7Ris(f )I (10) 
and Eoaalt es Fall 28E calf )+7E,s(f )I (1) 


in (10) is our proposed mixed tule. 
Theorem-3.2 
Let the smooth function f(p) is defined on —l<p <1, then the error cca s(f ) is associated 


with the mixed quadrature rule Reis (f is given by: 


f*(0) 


463 
150910, 


Ile 


Ear) 





Proof: 
The proof follows from (11). 
Theorem-3.3 


The truncated error bound for Ecous(f) =I(f)- Reons (f) is evaluated by: 


64M ix 
Beaualt \s 425258! where M =maxtf (p ) 
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Proof: 


i 
~ 12608! 


ae viii 
Exs (f ) = foal (7, ), hh € [- 11] 


Eca (f ) f "(n, ), mE [- 1,1] 


= 2 viii _ pviii 
Ecault) = gana lf (7) f (n,)| 


where K = max 


—-I<x<l 


[- 1, 1| such that § = fc). 





iy ma) and S = min 


—I<x<li 





Poe Hence there exists points d and C in the interval 


L 


ee 


4 
for some yel- 1, 1, 
where ld —¢ <2 


64 


<Dssal 


Ecens(f) 


64M ix 
Hence Focal } <79595<8) where M =max|f (p). 


4. NUMERICAL VERIFICATION 
In this section, Table 1 and Figure 2 for the numerical approximation of RLC-circuit for the mixed 


quadrature rule Roa 15 (f) with the initial conditions [ (0) =, 7 (0) = lare prepared. 


Table 1. Comparison of exact result with approximate result 


f(t) er ee lError| 

sint 0.841470984807897 amp 0.84 1492252269964 amp 0.00002 1267462067 

cost 0.41410934759 113 1 amp 0.41528766082 175 1 amp 0.001 178313230620 
e 1.041865355098910 amp 1.04183440 1921017 amp 0.00003095 3177893 


Comparison of Exact result 
with approximate 





Figure 2. Comparison of exact result with approximate 
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5. CONCLUSION 

The RLC circuit problem is resolved through integral equation of second kind with suitable kernel 
which is Volterra in nature. The approximate solution can be achieved by taking a few calculations which 
is much better. The results obtained by the presented method to introduce KVL for solving the said integral 
equation reveals found tobe more effective and convenient. The efficiency of the method is verified 
numerically by taking some test problems in Table 1 and it is depicted in Figure 2. The extensive work 


for finding instantaneous current for any value of R, L and C and for any function f(t) can be carried out. 
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